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Abstract —In this study, a comparison was made between 
experimental and numerical analysis results using the 
classical dam-break test case over dry bed. DualSPHysics 
software based on Smoothed Particle Hydrodynamics 
(SPH) method was used to make the numerical analysis. 
Experimental data were obtained from Kocaman [l]’s 
laboratory setup through image processing technique. It 
was observed that the numerical and experimental results 
are in good agreement. 
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I. INTRODUCTION 

Dams can be built to supply potable water, generate 
electricity and control flood. Since the massive volume of 
water may be found upstream, the dams pose risks for the 
residential areas located downstream. In case dam-break 
case, the massive volume of water can create a flood and 
this can cause many losses of life and property. The flood 
caused by the dam-break is more destructive than a rain 
induced flood. Agricultural and settlement areas are often 
those that have experienced floods caused by dam-breaks in 
the past. Investigating the dam-break problem scientifically 
is remarkable due to insufficient data of the eventuated 
dam-break cases. This and this kind of study aims at 
combining experimental data and numeric analysis to 
predict the areas that will be affected by the flood and 
contribute to the improvement of the early warning systems. 
Due to difficulties in obtaining real case-data, and although 
the hydraulic problems are large scale in the field, the data 
of the numerical analysis can be validated through small- 
scale experimental studies [2-5]. Most of the early 
numerical studies used Shallow Water Equations (SWE) to 
describe dam-break wave propagation [6-8]. Due to recent 
developments in computer power, it is possible to widely 
use numerical methods based on Computational Fluid 
Dynamics (CFD). Some researchers used SWE and 


Reynolds Averaged Navier-Stokes (RANS) equations to 
investigate dam-break flow in 2D and 3D [9,10]. 

Smoothed Particle Hydrodynamics (SPH) method is one of 
the most commonly used numerical methods in Hydraulics. 
The method has become popular recently. It was developed 
to investigate the astrophysics cases at the middle of the 
1970s [11]. 

SPH method is still developing and it is frequently used in 
the modelling of various engineering problems such as 
hydraulic [12,13], structural dynamics [14] and solid 
mechanics [15]. One of the common software used to model 
the problems is DualSPHysics. This freeware software is 
used to study free-surface flow phenomena where Eulerian 
methods can be difficult to apply, such as dam-breaks, 
sloshing, wave impacts on off-shore structures [16-17]. 

In this study, the numerical and experimental results were 
compared for the dam-break flow over a dry bed at initial 
stages. The comparison was made through water surface 
profile curves. The numerical results were obtained from 
the SPH method, whereas experimental results were 
determined using image processing technique. 

II. EXPERIMENTAL SET-UP AND NUMERICAL 

MODEL 

2.1. Experimental Set-up 

The experimental studies were carried out at Civil 
Engineering Hydraulics Laboratory of the Cukurova 
University, Turkey. The test setup was conducted by 
Kocaman [1]. Experimental studies were performed in a 
rectangular channel. Lateral and bottom surfaces are glass. 
Dimensions of the channel are 8.90 m x 0.3 m x 0.35 m 
(length, width, height, respectively). The plate representing 
the dam was placed at 4.65 m from the channel entrance and 
this part was considered as a reservoir. The water depth in 
the reservoir was taken 0.25 m. A red dye was added into 
the reservoir to better determine the free surface levels from 
recorded images. The downstream side of the channel is 
kept dry at the beginning of the experiment. The mechanism 
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illustrated in Figure 1 allowed for the instantaneous removal 
of the vertical plate to model the dam-break. The 4-mm 
thick plate was used as a dam. A steel rope was attached to 
the plate top; the rope was drawn over a pulley of 15 kg 
weight attached at the other end. By releasing the weight 
from 1.50 m above the floor, the plate was removed 
instantaneously. The plate removal time was determined as 
approximately 0.07 s from video records. 



Experimental results were obtained through image 
processing technique. Three adjacent cameras were 
simultaneously used to record the experiment. The recorded 
videos were calibrated to obtain the water surface profiles. 
The cameras are JVC c920e with acquisition rate 25 frames 
per second and resolution 752x586 pixels. A complete 
description of the experimental setup and measurement 
techniques can also be seen in Ozmen-Cagatay and 
Kocaman [9-18]. 

2.2. Numerical Model 

Dual SPHysics is a numerical model based on the SPH 
method. This software is capable to study real engineering 
problems by using graphical processing units (GPU) [16]. 
In this model, fluid, boundary and solid bodies are 
represented with particles. The interaction of the 
neighboring particles depends on the distance between the 
particles. A kernel function (W) defines to calculate the 
interaction, which is affected by smoothing length (h). A 
particle is involved in an interaction with the particles 
located in a circular area, which size depends on the 
smoothing length. There is no interaction between particles 
placed out of the circular area. Physical quantities such as 
pressure, location, density, velocity etc. are calculated 
through interpolation depending on the smoothing length. 
The numeric foundation of the SPH are based on the 


interpolation of integrals. An F function is defined as shown 
in Equation 1 [12]: 

F( r) = / F(r ') W (r - r', h)dr (1) 

where F(r) is the average, W is the kernel function, h is the 
smoothing length. If this function is interpolated in an “a” 
particle, Equation 2 is obtained; 

F(r a ) * F(r b ) W(r a - r b , h)Av b (2) 

where, subscripts a and b are individual particles. Av b is the 
volume of the neighboring particle b. Since the volume is 
the ratio of the density to mass Equation 2 becomes [19] 

F (r a ) * ^ F(r ^ © W(Ta ~ Tb ' h) (3) 

where, nib and pb are the mass and the density of the particle 
b respectively. 

The accuracy of an SPH analysis depends on the kernel 
function. The function is defined as q=r/h. Here, r is the 
distance between two particles and q is the dimensionless 
ratio. The definition of the Wendland kernel used in this 
study is given in Equation 4: 

W(r, h) =a D (l- 0 (2 q + 1) (4) 

where old value is 7/4nh 2 for 2D analysis while, 21/1 6nh 3 for 
the 3D analysis. 

In the SPH formulation, the Navier-Stokes equations are 
solved based on the assumption of weakly compressible 
fluid [20]. Momentum equation is used to define the 
acceleration for an “a” particle with interaction of the 
neighboring particles: 

where, v is velocity, P is pressure, g is the gravity and W a b is 
kernel function. TJ a b is a term based on the artificial 
viscosity created by Monaghan [20] 

Since the mass of the particles is constant, the conservation 
of mass or continuity equation is used to measure the 
density variations: 

‘% L = Z. b m b v ab .V a W ab (6) 

The relationship between pressure and density is given by 
the Tait’s equation of state, shown in Equation 7 [21]: 

P = (7) 

b = c$p 0 /Y ( 8 ) 

where y = 7 and po= 1000 kg/m 3 . 
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The speed of sound c 0 is numerically computed like at least 
ten times the maximum velocity in the system. For the dam- 
break flow the wave-front velocity was used as follows: 

Co = CsyfgK, (9) 

where c s is the coefficient sound, h w is the still water depth. 
If the value of c* is increased, compressibility of the fluid is 
decreased according to Equations 7-8. 

The model parameters set in the DualSPHysics software for 
the SPH method to make the numerical analysis, are 
reported in Table 1: 


Table. 1: Parameters of the numeric model 


Artificial viscosity 

0.01 

Courant coefficient 

0.2 

Coefficient Sound 

20 

Distance of the particles 

0.002m 



DeltaSPH value 

0.1 

Number of the particles 

173650 

Smoothing length 

0.002828 

Step Algorithm 

Symplectic 

Kernel 

Wendland 


III. VALIDATION OF THE NUMERICAL MODEL 

Comparisons between the SPH numerical and experimental 
free surface profiles at initial stages can be seen in Figure 2. 
It is observed that when the plate moves above, the water 
volume in the upstream rapidly propagates towards 
downstream. Parabolic wave front at the initial stage 
becomes hump shaped as time progressed with surface 
roughness. It can be understood that the surface roughness 
affects the results at the lately stages more than initial 
stages. The SPH results agree well with experiment. 


■ 




(a) (b) 

Fig. 2: Comparison of free surface profiles by: a) SPH b) experiment 


Figure 3 shows the comparison between the computed free 
surface profiles using the SPH and the experimental data at 
different times. The flow depths h and horizontal distances v 
were non-dimensionalized by the initial water depth ho. 
Similarly, time t was multiplied by (g/ho) 1/2 to obtain 
dimensionless time. A good agreement is observed between 
the measured and SPH-computed free surface profiles at all 
times. However, the wave-front velocities determined by 
the SPH model are always higher than experiment except at 
very early stage (T=1.13). 

It was observed that slight changing of some parameters 
(number of the particles, distance of the particles, artificial 
viscosity, smoothing length etc.) have remarkable effects on 
the numerical results. Especially, the particle resolution in 
the SPH has a great impact on the accuracy and stability. 


Reducing the distance of the particles causes an increase in 
the number of the particles. For this reason, solving abilities 
of the computer becomes insufficient. High-performance 
computing (HPC) needs to be analyzed for the numerical 
model with a large amount of particles. However, GPU 
accelerated computing gives advantages to the SPH method 
compared to other CFD based numerical methods. 
DualSPHysics has the GPU option which provides a higher 
computing power than CPU in order to accelerate SPH 
simulation with low cost. In this study, runtime of the 
problem took 556 s for 1.5 simulation time by using 
Geforce GTX 580 card with 3GB RAM. It was also 
considered that the variations of the parameters should be 
investigated in further studies by using different 
experimental data having obstacles, wave impacts, etc. 


www.ijaers.com 


Page | 211 




















International Journal of Advanced Engineering Research and Science (IJAERS) 
https://dx.doi.Org/10.22161/iiaers.4.12.31 


[Vol-4, Issue-12, Dec- 2017] 
ISSN: 2349-6495(9) / 2456-1908(0) 







Fig. 3: Experimental and numeric results of water surface profiles 


IV. CONCLUSIONS 

In this study, the flood wave resulted from dam-break 
phenomenon over a dry bed at initial stages was examined 
by comparing the SPH based numerical model with 
experimental one. A very good agreement was observed for 
free surface profiles. However, it was found that wave-front 


velocities of the SPH model is higher than the experimental 
one. 

Since the results are in good agreement, it can be concluded 
that DualSPHysics based on the SPH method can be 
utilizable and useable to analyze the chaotic hydraulic 
problem such as dam-break. However, in order to be able to 
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apply large-scale real-life problems, further validation is 
required with experimental studies representing different 
topography conditions. 
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